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A magnetic levitation system (MLS) is a complex nonlinear system that 
requires an electromagnetic force to levitate an object in the air. The 
electromagnetic field is extremely sensitive to noise which can cause the 
acceleration on the spherical object, leading it to move into the unbalanced 
region. This paper presents a comparative assessment of controllers for the 
magnetic levitation system using proportional integral derivative (PID) 
controller based optimal tuning. The analysis was started by deriving the 
mathematical model followed by the implementation of radial basis function 
neural network (RBFNN) based metamodel. The optimal tuning of the PID 
controller has offered better transient responses with the improvement of 
overshoot and the rise time as compared to the standard optimization 
methods. It is more robust and tolerant as compared to gradient descent 
method. The simulation output using the radial basis based metamodel 
approach showed an overshoot of 9.34% and rise time of 9.84 ms, which are 
better than the gradient descent (GD) and conventional PID methods. For the 
verification purpose, a Simscape model has been developed which mimic the 
real model. It was found that the model has produced about similar 
performance as what has been obtained from the Matlab simulation. 
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1. INTRODUCTION 


For decades a % was considered as a dream in science fiction, in 1726, Jonathan Swift described the 
magnetic levitation system (MLS) for the first time. Also, in 1842 an English clergyman called Samuel 
Earnshaw described the importance of maglev and its limitations. He showed that the system of maglev has 
instability issues, where the force between the static magnets and the contactless levitated part was 
impossible to be stable. The free levitated part has unstable displacement in at least one direction [1], [2]. 
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Maglev has many applications in engineering field and industrial sector, such as in magnetic bearing 
systems, high-speed trains, vibration insulation systems, stepper photolithography, and wind turbine [3]. 
Recently, magnetic levitation systems have been appreciated for removing mechanical contact friction, to 
reduce maintenance costs and to achieve high-precision positioning. In 2010, a group of researchers from 
Rice University had developed a three-dimensional tumor model related to magnetic levitation. They have 
injected the cancer cells with magnetic iron oxide and gold nanoparticles. Moreover, by installing a coin size 
magnet near the infected area, they have successfully lifted the cells [4]. 

Electrodynamic suspension (EDS) and, electromagnetic suspension (EMS) are the two types of 
magnet levitation methods, as illustrated in Figure 1. “Repulsive levitation” is a term used to describe EDS 
systems. To supply comparable levitation sources, superconductivity magnets [5] or permanent magnets [6] 
are utilized. The magnetic levitation force provided by EDS is relatively steady and allows for a large 
clearance. However, as compared to the EMS method, the magnetic materials production process is more 
difficult and costlier. The EMS system referred to attractive levitation system. Since magnetic levitation force 
is inherently unstable, regulating the system is considerably more difficult than controlling the EDS system, 
and extra power is needed to maintain levitation height [7], [8]. 
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Figure 1. Attractive vs repulsive maglev system 


Over the years, engineers and researchers have been paying great attention to stabilize the magnetic 
levitation system. The characteristic of the maglev system is extremely nonlinear. It is unstable and 
considerable as uncertainty. Proportional integral and derivative (PID) has been utilized to control the maglev 
system, which for so many years has been used in the industrial field. PID controllers have been utilized 
since the 1890s for controller design [9]-[12]. Until today the industrial field still uses PID controllers with 
other optimization techniques. The PID controllers could stabilize the maglev system, although the control 
performance of the system is limited due to the fixed controller parameters [13]-[16]. There are many optimal 
techniques that have been used to optimize PID parameter and one of these is neural network, and the 
increase of neurons may affect simulated response, because of the computation process. The metamodeling 
approach which will be explain in section 3, intended to reduce the computation process and time process, 
besides produce more efficint response. The research method will be presented in section 2, followed 
controller design in section 3, result, and discussion in section 4, and finally the paper conclusion in 
section 5. 


2. METHOD 

Figure 2 shows the basic construction of the maglev system. Where the ferromagnetic ball is 
levitated in a controlled electromagnetic field [13]. The system comprises of a position sensor, 
electromagnetic actuator, optimal PID controller and a power supply circuit. The system’s parameters used 
for this study are listed in Table 1. By studying the mechanical and mathematical modelling, behavior of the 
system can be modelled as in the following subsections. 
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Figure 2. Magnetic levitation model 


Table 1. Real parameters of maglev system 


Parameters Values Units 
M 0.5 Kg 
R 20 Ohm (Q) 
L 0.158 H 
g 9.8 m/sec? 
K 0.08 
Y 0.03 m 
i 1.44 A 
K, -458.7157 
K, 5.8929 


2.1. Electrical equations of the MLS 
Magnetic levitation system consists of resistor and inductance, and with the current flow through the 


circuit, it produces a magnetic field, the energy needed to levitate the Ferro ball, thus, to calculate the energy 
stored inside the magnetic field will be as: 


p= iv (1) 
di 
v=1t 2) 


the current starts at zero to steady-state at I final, will get: 
U= fav = Lf! idi 
1 
t=, Li? (3) 


2.2. Mechanical equations of the MLS 
By utilizing the Newton’s third law of motion, and by ignoring the damping force and friction force 


of air [14], [15], the formula will be as (4) and (5): 
mx = mg - f (Xp, ic) (4) 
ss ic? 
mx = mg- kee (5) 
where: m is the mass of the ball, k is magnetic force constant, x is ball position and g is gravity constant. 


2.3. Transfer function model of the MLS 
The electromagnetic force F(i,x) is nonlinear, and thus the MLS is a typical nonlinear system 


which requires linearization [16], as shown by (4) and (5). Linearization is achieved by utilizing Taylor's 
method to extend the nonlinear function. At the equilibrium point of ig and x9, Taylor's theorem is applied to 
the maglev system in (5) will yield, 


F(i, x) = F(io, Xo) + rG (i- io) + PAL) g- xo) (6) 
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when the ball in maglev system in equilibrium condition the (4) and (5) will become: 
mg + Fhios Xo) =0 & mg = —F(ig,X>) (7) 


by defining the maglev system input will be the voltage to the inductor (Ka), and the feedback from the plant 
will be the sensor gain (K,). By substituting the parameter’s values from Table 1 will yield: 


77.84 


GS) = cous? 2083 A 
Converting the transfer function in (8) to state space form will yield. 
0 1 0 
~ peas 0 Bee aket 1], D= [0] (9) 


2.4. PID controller based on ziegler nichols tuning 
PID known as "proportional-integral-derivative controller" is a linear controller and a typical 
method used in industrial application [17], [18]. PID controller is mathematically represented by: 


u(t) = K,e(t) +K; f e(t)dt + Kp ae (10) 
where e(t) is the error. Ziegler Nichols based PID [19]. Equation is given as (11) and (12): 
u(t) = Ky (ece) +z f e(a)dt+ Tq w) 
= K,e(t) +2 fy e(t)dt + Kp Ta ŽP (11) 
Ka = KpTa. K; = To Kp = 0.6 X Ky (12) 


where, K, is the ultimate gain and Tg and T; is the derivate time and integral time correspondingly. 


2.5. Radial basis function (RBF) with gradient descent PID tuning 

In the optimization studies, there is no definitive answer as to which model is better than the others. 
However, the Kriging model, Gaussian, and radial basis function are still the most often used methods [20]. 
The general radial basis function (RBF) equation is given as in (13): 


Y(k) = X well — Gill) U2) 
where, x e R®*+ is an input vector, ¢ is a basis function, C; is the RBF centres in the input vector space, N is 
the number of hidden layer neurons, and w; is the weights in the output layer [21]. 

Gaussian function usually called Gaussian and used to arbitrary real constant and non-zero variable; 
it is simple and using it with the RBF network does not increase the computational complexity when using it 
with a significant variable or large dataset. where better analysis results can be obtained when executing the 
data set [22], [23]. The Gaussian function expressed as (14): 


af 
h; =exp exp E ei |i =1,2,..,m (14) 


2b} 
the weight of each vector is found as (15): 
W = [w;]",j = 12,3... N (15) 
the identified RBF network output is: 
Ym (k) = Diy Wj (k)h;(k) (16) 


the identifier for performance cost function is: 


An optimal proportional integral derivative tuning for a magnetic levitation ... (Abdualrhman Abdalhadi) 


1360 O ISSN: 2502-4752 
IE) =} (YE) = ym)? (17) 


According to the network model chosen above, the iterative method is used to determine the updated 
weights, w, biases, b and control, c. to be calculated as: 


w,(k) = wlk — 1) + n(Y) — ym) aj) + a (w;(k -1)- wk - 2)) (18) 
AbjCK) = (Y) = yn OE) we) BoE 
bj(k) = by (Ck 1) + nAby(k)) a (bj(k - 1) - yk - 2)) (19) 


Abjilk) = (Y) — Yn (I))w; (ke) FE- 


bj(k)? 


cll) = cp (R= 1) + Ac (Ie) a (jek - 1) — cilk- 2)) (20) 


where: 77 is the learning rate and æ is the momentum factor. 

The Jacobian matrix is used to compare output values and evaluate the system's control sensitivity 
by executing the output and input of the controlled system in real-time [24], [25]. The Jacobian algorithm is 
provided as: 


dy(k) _ A¥m(K) _ vm cjil(k) 
aaua © aaucey — 2-1 iFM) oe (21) 
based on the algorithm of incremental PID to the control system is given as: 
e(k) = r(k) — y(k) (22) 
the inputs of the PID controller are given as: 
xe, = e(k) —e(k — 1) (23) 
xc, = e(k) (24) 
xc3 = e(k) — 2e(k — 1) + e(k — 2) (25) 
the control algorithm is taken from PID output as: 
Au(k) = kpxcı + kixc, + kgxcs (26) 
three parameters of PID controller used to adjust the gradient method: 
ðE a 
AK, = -135 = ne(k) xc, (27) 
ðE a 
AK; = —n 5 = nelk) = xc (28) 
ðE a 
AK4 = =n 5-7 = ne(k) 5 xcs (29) 


the PID parameter increment AKp, Aki, AKd can be updated through the iteration. The parameters of PID 
controller kp, ki, ką can be written as: 


kp(k) = kp(k — 1) + Akp (30) 
ka(k) = ka(k — 1) + Aka (32) 
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2.6. RBFNN metamodel based PID tuning 
Metamodelling process involves the determination of simpler model from a complex model that 

provides less computation or control effort. Table 2 shows the RBF metamodel parameter for both the initial 

dataset and the large dataset. The metamodel approach used with radial basis function neural network is 

summarized in the following steps: 

— Select the parameters need to be optimized (i.e., PID controller parameters). 

— Next, generate initial input dataset (D) with a smaller number of iterations for the PID controller 
parameters. 

— Run the simulation in simulink (as shown in Figure 3) by using all input dataset (D) and set the target 
dataset (T) from the corresponding integral square error (ISE) index values. 

— Run the fitting of D and T dataset using radial basis function neural network (RBFNN). 

— Evaluate RBFNN obtained in (4) using denser dataset (large dataset, D’). 

— Find the minimum output of the RBFNN output in (5) which corresponds to the smallest value of the ISE 
for large dataset. 

— The equivalent PID gains that minimized the RBF output, will be the gains that will be used to simulate 
the maglev system, which is finally identified as the PID optimal parameters. 

— The integral square error (ISE) block is connected to magnetic levitation system MLS shown in Figure 3 
used to check the error, thus the process will be repeated if the output min error is not satisfied. In other 
word, the error gotten by metamodel RBF must be less than the one gotten by conventional PID. 


Table 2. Controller parameter used for simulation 


Dataset PID Parameter PID Gain Range 
Large data set Ki 4: 0.1: 20 
K, 50: 2: 120 
Kp 0.06: 0.01: 0.40 
Total data configuration 202,860 
Initial data set Ki 9: 0.1: 15 
K, 50: 2: 116 
Kp 0.06: 0.01: 0.13 
Total data configuration 16,592 
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Figure 3. Overview of metamodel based RBF 


3. RESULTS AND DISCUSSION 
3.1. RBFNN metamodel-based PID tuning performance 

The metamodeling approch intended to reduce the computation process and time process, besides 
produce more efficint response Figure 4 shows the magnetic levitation model implementation in 
Matlab/Simulink. The varying parameters of PID values (i.e., represents as inputs in RBFNN metamodel) in 
different sampling interval were investigated to see the resultant series of ISE values. 

The RBFNN metamodel based PID tuning was used to get the best performance of maglev system 
by obtaining the suitable PID gains. The system response is shown in Figure 5. where the PID values used 
are: kp = 15, k; = 100,kq = 0.1 when the ISE during the tuning was 0.00104. The system output 
characteristics are, rise time T, = 9.85 ms, overshoot = 9.34%, undershoot = 0.094%, settling time T; = 
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0.25 s. comparing these results with the one in [3], [14]. The rise time is 0.28 sec for two degrees of freedom 
controller and the settling time is 0.35sec. also comparing 0.25 sec settling time to [16], the best settling time 
they got is 0.6 sec. the results for coefficient diagram method (CDM) based PID controller in [18] shows a 
good overshoot, but slow settling time of 0.96 sec. Wahid et al. [24] RBF based PID controller has been used 
for maglev system and the results show a settling time of 0.5 sec. 
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Figure 5. Response of RBFNN metamodel based PID tuning 


3.1. Gradient descent based PID tuning performance 

The gradient descent (GD) is an algorithm of optimization which is intuitive and efficacious. In 
many diverse electromechanical systems, the GD technique has effectively optimized the complicated 
system. The method used in this study to tune PID controller is the gradient descent provided by 
Matlab/Simulink toolbox, which called as response optimizer. The output response is shown in Figure 6. The 
PID parameters yield from GD optimizer are kp = 7.696, k; = 51.34, kg = 0.46. The system output after 
applying the PID gain it show that, rise time T, = 72.92 ms, overshoot = 21.34%, undershoot = 1.91%, 
settling time T, = 0.70 s. 
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Figure 6. Response of GD Based PID Tuning 


3.2. Comparative study of tuning performance 


2.5 


The maglev system response as shown in Figure 7 has three different tuning methods which are the 
RBFNN metamodel, gradient descent and conventional PID. Besides, the comparison data between these 
tuning methods are listed in Table 3. The best rise time is 9.85 ms and the overshoot is 9.34%, which is 
offered by the RBFNN metamodel approach. However, the settling time performance is very close between 
conventional PID and RBF metamodel. The conventional PID has better results than the gradient descent 
method where the rise time and the overshoot are 18.45 ms and 19.53%, respectively, thus the less successful 
method used to tune the PID controller in this study has been found to be produced by the GD algorithm. 
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Figure 7. PID, meta-RBF and GD output response 


Table 3. Performance comparison of tuning methods 
Response characteristic | Metamodel GD Conventional tuning 


Rise time (ms) 9.85 72.92 19.535 
Overshoot (%) 9.34 21.34 18.45 
Undershoot (%) 0.144 1.919 1.553 
Settling time (s) 0.25 0.70 0.22 


2.5 
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3.3. Performance verification of maglev system in Simscape simulation 

Simscape multibody is a library integrated with Matlab/Simulink and uses blocks to build a 3D 
model that consists of joints, constraints, and force elements for instance electromagnetic force. The 
developed block connection of Simscape for maglev system is shown in Figure 8. After connecting all the 
blocks together and setting the parameters value based on the actual model, the maglev system is then linked 
to 3D animation design via solid works. There are two important components for Simscape to work, first, is a 
compatible software that can design 3D model and export that model to a Simscape file, and here it comes 
the role of solid works, where solid works computer-aided design (CAD) can design separate parts based on 
specific dimensions and then sample them all in one 3D model to be exported as a Simscape file. The second 
part is reading the Simscape file using Matlab/Simulink and define all the parts based on their functionality. 
The 3D design is shown in Figure 9. 
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Figure 9. Maglev 3D diagram using solidworks 


The values for PID parameters that were used to verify the maglev system using Simscape model 
are kp = 15,k; = 100 and kg = 0.1, which utilizes RBFNN metamodel optimal parameters, as obtained in 
section 4. The system output performance is depicted in Figure 10. The system produces an overshot of 
13.07% and a rise time of 7.38 ms. It is found that the performance of the simulated system (as in section 4) 
has better results in terms of overshoot (i.e., 9.34%), as the previous simulation represented the ideal 
environment. However, for the Simscape model, we build a physical component based on physical 
connection which is representing a real system, hence, the system performance by using similar PID 
parameters yield a slightly reduced performance, especially in term of the overshoot. 
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Figure 10. Output response for MLS Simscape model 


4. CONCLUSION 

This paper praises a PID tuning method using the radial basis function neural network (RBFNN) 
Metamodel and comparing the output results with gradient descent optimization and conventional PID 
methods. Furthermore, the RBF metamodel controller has been compaered with different papers that can be 
seen in the results section, and the compareson shows that the method used in this paper has improved the 
performance of settling time and overshoot between 35% to 70%. The metamodel based PID tuning approach 
showed transient response with better overshoot and rise time as compared to other optimization methods. 
For the trained datasets, RBF networks perform more robustly and indulgently than the GD method even 
when dealing with noisy input data sets. Thus, RBFNN metamodel is recommended specifically for surface 
with regular peaks and valleys for functional approximation issues, since effective and accurate designs are 
attainable. The verification using Simscape model had also shown that the tuned PID parameters using the 
proposed approach able to obtain about similar performance as of using Matlab simulation. 
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